Tarea

Emily B - C01213

Carga de paquetes

library(sf)
library(raster)
library(dplyr)
library(spData)
library(leaflet)
library(plotly)
library(DT)

Carga de datos

primates <-
  st_read("https://raw.githubusercontent.com/gf0604-procesamientodatosgeograficos/2021i-datos/main/gbif/primates-cr-registros.csv",
    options = c(
     "X_POSSIBLE_NAMES=decimalLongitude",
      "Y_POSSIBLE_NAMES=decimalLatitude"),
    quiet = TRUE
  )

# CRS
st_crs(primates) = 4326

# Capas de cantones y distritos
cantones <-
  st_read(
    "https://raw.githubusercontent.com/gf0604-procesamientodatosgeograficos/2021i-datos/main/ign/delimitacion-territorial-administrativa/cr_cantones_simp_wgs84.geojson",
    quiet = TRUE
  )
provincias <-
  st_read(
    "https://raw.githubusercontent.com/gf0604-procesamientodatosgeograficos/2021i-datos/main/ign/delimitacion-territorial-administrativa/cr_provincias_simp_wgs84.geojson",
    quiet = TRUE
  )

Obtencion del nombre del canton

primates <-
  primates %>%
  st_join(cantones["canton"])

Introduccion

A continuacion se presentan una serie de datos estadisticos mediantes graficos y un mapa, de cuatro diferentes especies de primates que habitan Costa Rica, esto con la finalidad de mostrar y visulizar su ocupacion del pais.

Tablas de registros de presencia

Grafico de pastel

Filtros

mono_arana <- primates %>%
  dplyr::select(species,
         stateProvince,
         locality,
         eventDate,
         decimalLatitude,
         decimalLongitude) %>%
  filter(species == 'Ateles geoffroyi')

mono_carablanca <- primates %>%
  dplyr::select(species,
         stateProvince,
         locality,
         eventDate,
         decimalLatitude,
         decimalLongitude) %>%
  filter(species == 'Cebus capucinus')

mono_aullador <- primates %>%
  dplyr::select(species,
         stateProvince,
         locality,
         eventDate,
         decimalLatitude,
         decimalLongitude) %>%
  filter(species == 'Alouatta palliata')

mono_ardilla <- primates %>%
  dplyr::select(species,
         stateProvince,
         locality,
         eventDate,
         decimalLatitude,
         decimalLongitude) %>%
  filter(species == 'Saimiri oerstedii')

Capa de altitud

altitud <- getData(
  "worldclim",
  var = "alt",
  res = .5,
  lon = -84,
  lat = 10
)
altitud <-
  altitud %>%
  crop(provincias) %>%
  mask(provincias)

Colores de altitudes

 colores <- colorNumeric(
    c("#99FF66", "#FFFF09", "#E02020"), 
    values(altitud), 
    na.color = "transparent"
  )

Mapa

arana <- paste0("<b>","Especie: ","</b>",
      (mono_arana$species), "<br>",
      "<b>","Provincia: ","</b>",
      (mono_arana$stateProvince), "<br>",
      "<b>","Canton: ","</b>",
      (mono_arana$locality), "<br>",
      "<b>","Fecha: ","</b>",
      (mono_arana$eventDate), "<br>",
      "<b>","Coordenadas: " ,"</b>",
      (mono_arana$decimalLongitude), "<br>",
      (mono_arana$decimalLatitude)
      )
ardilla <-paste0(
  "<b>","Especie: ","</b>",
      (mono_ardilla$species), "<br>",
  "<b>","Provincia: ","</b>",
      (mono_ardilla$stateProvince), "<br>",
  "<b>","Canton: ","</b>",
      (mono_ardilla$locality), "<br>",
  "<b>","Fecha: ","</b>",
      (mono_ardilla$eventDate), "<br>",
  "<b>","Coordenadas: " ,"</b>",
      (mono_ardilla$decimalLongitude), "<br>",
      (mono_ardilla$decimalLatitude)
      )
aullador <-paste0(
  "<b>","Especie: ","</b>",
      (mono_aullador$species), "<br>",
  "<b>","Provincia: ","</b>",
      (mono_aullador$stateProvince), "<br>",
  "<b>","Canton: ","</b>",
      (mono_aullador$locality), "<br>",
  "<b>","Fecha: ","</b>",
      (mono_aullador$eventDate), "<br>",
  "<b>","Coordenadas: " ,"</b>",
      (mono_aullador$decimalLongitude), "<br>",
      (mono_aullador$decimalLatitude)
      )
carablanca <-paste0(
  "<b>","Especie: ","</b>",
      (mono_carablanca$species), "<br>",
  "<b>","Provincia: ","</b>",
      (mono_carablanca$stateProvince), "<br>",
  "<b>","Canton: ","</b>",
      (mono_carablanca$locality), "<br>",
  "<b>","Fecha: ","</b>",
      (mono_carablanca$eventDate), "<br>",
  "<b>","Coordenadas: " ,"</b>",
      (mono_carablanca$decimalLongitude), "<br>",
      (mono_carablanca$decimalLatitude)
      )

Mapa

primates %>%
leaflet() %>%
  addProviderTiles(providers$OpenStreetMap.Mapnik, group = "OpenStreetMap") %>%
  addProviderTiles(providers$Stamen.TonerLite, group = "Stamen TonerLite") %>%
  addProviderTiles(providers$Esri.WorldImagery, group = "Imagenes ESRI") %>%
  addRasterImage(
    altitud, 
    colors =  colores, 
    opacity = 0.5,
    group = "Altitud"
  ) %>%
 addCircleMarkers(
    data = mono_arana,
    stroke = F,
    radius = 2,
    fillColor = '#056718',
    fillOpacity = 1,
    popup = arana,
    group = "Ateles geoffroyi"
    )%>%
     addCircleMarkers(
    data = mono_ardilla,
    stroke = F,
    radius = 2,
    fillColor = '#35427B',
    fillOpacity = 1,
    popup = ardilla,
    group = "Saimiri oerstedii"
    ) %>%
     addCircleMarkers(
    data = mono_aullador,
    stroke = F,
    radius = 2,
    fillColor = '#A40C0C',
    fillOpacity = 1,
    popup = aullador,
    group = "Alouatta palliata"
    ) %>%
     addCircleMarkers(
    data = mono_carablanca,
    stroke = F,
    radius = 2,
    fillColor = '#B98E4F',
    fillOpacity = 1,
    popup = carablanca,
    group = "Cebus capucinus"
    ) %>%
  addLayersControl(
    baseGroups = c("OpenStreetMap", "Stamen TonerLite", "Imagenes ESRI"),
    overlayGroups = c("Ateles geoffroyi","Cebus capucinus","Alouatta palliata","Saimiri oerstedii", "altitud")
  ) %>%
      addMiniMap(
    tiles = providers$Stamen.OpenStreetMap.Mapnik,
    position = "bottomright",
    toggleDisplay = TRUE
  )